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In the Bak-Sneppen (BS) evolution model [1], random 
numbers, /j, chosen from a flat distribution between 
and 1, p(f), are assigned independently to each species 
located on a c?-dimensional lattice of linear size L. At 
each time step, the extremal site, i.e., the species with 
the smallest random number, and its 2d nearest neigh- 
boring sites, are assigned 2d + 1 new random numbers 
also chosen from p(f). This updating continues indefi- 
nitely. After a long transient process the system reaches 
a statistically stationary state where the density of ran- 
dom numbers in the system is uniform above f c (the self- 
organized threshold) and vanishes for / < f c . 

Despite the fact that it is an oversimplification of real 
biological process, BS model exhibits such common inter- 
esting features observed by paleontologists [2,3] as punc- 
tuated equilibria, power-law probability distributions of 
lifetimes of species and of the sizes of extinction events. 
These behaviors suggest that the ecology of interacting 
species might have evolved to a self-organized critical 
state. 

BS model displays spatial-temporal complexity, which 
also emerges from many natural phenomena, such as frac- 
tals [4], 1// noise [5], etc. This strongly suggests that 
various complex behaviors may be attributed to a com- 
mon underlying mechanism. Authors in Ref. [6] suggest 
that the relation of these different phenomena can be 
established on the basis of their unique models. It is 
even proposed by them that spatial-temporal complexity 
comes out as the direct results of avalanche dynamics in 
driven systems, and different complex phenomena are re- 
lated via scaling relations to the fractal properties of the 
avalanches. It hence can be inferred that avalanche dy- 
namics plays a key role in dealing with complex systems, 
especially when one needs to know the macroscopic fea- 
tures of the systems, since lingering on the inner structure 
of individuals will not be helpful [7]. 

Avalanche is a kind of macroscopic phenomenon driven 
by local interactions. The size of an avalanche may be 
extremely sensitive to the initial configuration of the sys- 
tem, while the distribution of the sizes (spatial and tem- 
poral) of avalanches, i.e. the "fingerprint" should be ro- 
bust with respect to the modifications, due to the univer- 
sality of complexity and the definition of self-organized 
criticality (SOC) [8]. In this sense, the extent that we 
know about avalanche will determine to what extent we 
do a complex system. Avalanche dynamics provides in- 
sight into complexity and enables one to further investi- 
gate the system studied. 

Though avalanche dynamics may be a possible un- 
derlying mechanism of complexity, the definitions of 
avalanches can be vastly different for various complex 
systems, or for same sorts of systems, even for the same 
one. In BTW model [9], an avalanche is intrigued by the 
adding of a grain or several grains of sand into the sys- 
tem. The avalanche is considered over when the heights 
of all the sites are less than the critical value, say, 4. 
In BS model [1,6], several types of avalanches, for in- 
stance, /o-avalanche, G(s)-avalanche, forward avalanche 



and backward avalanche, etc, are presented. These dif- 
ferent definitions of avalanches may show their unique 
hierarchal structures, while they manifest the common 
fractal feature of the complex system, that is, SOC. It can 
be inferred that various types of avalanches are equiva- 
lent in the sense that they imply complexity. 

Since similar structures and common features evidently 
arise in different types of avalanches, it is straightforward 
that various avalanches differ each other only in the con- 
texts from which one comprehends them. As known, the 
major aim of avalanche study is to investigate the uni- 
versal rules possibly hidden behind the evolution of the 
systems or the models. Hence, the means of understand- 
ing the avalanches appear crucial. Better ways may en- 
able one to know more about the system or the model 
and hence to have better comprehension of the features 
corresponding to complexity. From this point of view, 
when studying avalanches one should try to choose the 
easier ways instead of the more difficult ones. 

The evolution of the highly sophiscated BS model 
shows a hierarchal structure specified by avalanches, 
which correspond to sequential mutations below certain 
threshold. It has been noted [10] that in BS model an 
avalanche is initiated when the fitness of the globally ex- 
tremal site (the species with the least random number) 
is larger than the self-organized threshold. That is, the 
triggering event of an avalanche is directly related to the 
fitness, the feature of individuals. In other words, the 
avalanche is directly associated with the feature of in- 
dividuals instead of general features of the ecosystem 
whole. Is it feasible that the avalanches are di- 
rectly intrigued by the global feature of the whole sys- 
tem? Can such global feature be expressed in terms of 
the corresponding quantity? If such quantity found and 
such avalanches observed, may the new avalanches pro- 
vide a new and easier way in investigating properties of 
the model? 

One of our previous works [10] presents such a different 
hierarchy of avalanches (/o-avalanche) for BS model. We 
defined a global quantity, /, which denotes the average 
fitness of the system. The new type of avalanches are 
directly related to /. In this paper, we present a mas- 
ter equation for the hierarchal structure of / - avalanches. 
It prescribes the cascade process of smaller avalanches 
merging into bigger avalanches when the critical param- 
eter /o is changed. An infinite series of exact equations 
can be derived from this master equation . The first or- 
der exact equation, together with an scaling ansatz of the 
average sizes of avalanches, shows the exact result of 7, 
the critical exponent governing self-organization, to be 
universally 1 for all dimensional BS models, which has 
been confirmed by extensive simulations of the model. 
We also establish scaling relations related to some criti- 
cal exponents for /o-avalanche and make predictions on 
the values of some exponents. 

The quantity / is a global one of the ecosystem and can 
be expected to involve some general information about 
the whole system. It may represent the average popu- 
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lation or living capability of the whole species system. 
Larger / shows that the average population is immense 
or the average living capability is great, and vice versa. 
/ is defined as 

L d 
i=l 

where fi is the fitness of the ith species of a system con- 
sists of L d species. Let BS model start to evolve. At 
each time step of the evolution, apart from the random 
numbers of the globally extremal site and its 2d nearest 
neighboring sites, the signal / is also tracked. Initially, 
/ tends to increase step-wisely. As the evolution contin- 
ues further, / approaches a critical value f c and remain 
statistically stable around f c . The plot of / versus time 
step s shows that the increasing signals of / follow a 
Devil's staircase [8] , which implies that punctuated equi- 
librium emerges. Denote F(s) the gap of the punctuated 
equilibrium. Actually, F(s) tracks the peaks in /. After 
some careful derivation one can write down an exact gap 
equation [6,10] 

dF(s) f c - Fjs) 
ds Li{S) F(s) ' 1 > 

where (S)p^ denotes the average size of avalanches oc- 
curred during the gap F(s) when / < F(s). This exact 
gap equation will be exactly solved in this paper. 

Signals f(s) play important roles in defining /q_- 
avalanche. For any value of the auxiliary parameter fo 
(0.5 < /o < 1.0), an /o-avalanchc of size S is defined as 
a sequence of S — 1 successive events when f(s) < fo 
confined between two events when f(s) > fo- This 
definition ensures that the mutation events during an 
avalanche are spatially and temporally correlated. It can 
also guarantee the hierarchal structure of the avalanches: 
larger avalanches consists of smaller ones. As fo is raised, 
smaller avalanches gather together and form bigger ones. 
The statistics of fo will inevitably have a cutoff if fo is not 
chosen to be f c . This will not affect the size distribution 
provided that fo approaches f c . Extensive simulations 
show that exponents r of /o-avalanche size distribution 
are 1.800 and 1.725 for ID and 2D BS models respec- 
tively, amazingly different from the counterparts of f - 
avalanche, 1.07 and 1.245 [6]. This strengthens the spec- 
ulation that /o-avalanche is a different type of avalanche, 
distinguished from any types of avalanches found previ- 
ously. 

Denote P(S, fo) the probability of acquiring a Jo- 
avalanche of size S. The signals f(s) (fo < f(s) < 
fo + dfo) will stop the /o-avalanches and not (fo + dfo)- 
avalanches. That is , as fo is raised by an infinitesimal 
amount dfo some of /o-avalanches merge together to form 
bigger (fo + d/o)-avalanches. This exhibits a hierarchal 
structure of /o-avalanches and will be prescribed by the 
below exact master equation. In some sense, the master 



equation reflects the "flow" of probability of avalanche 
size distribution with respect to the change in fo- 

Simulations show that / approaches f c and remain 
statistically stable in the critical state. This feature is 
greatly different from the feature of / m ;„ (fitness of glob- 
ally extremal site), which can vary between and 1. 
While / in the critical state fluctuates slightly around f c . 
Therefore, the / -avalanches will have no good statistics 
if fo is chosen as the value far less than f c , since there 
only exists smaller avalanches in the model. To acquire a 
better and reasonable distribution of / -avalanches sizes, 
one should choose the value of fo under the condition 
fo — > f c - It should be emphasized that the master equa- 
tion listed below is valid also for fo — > f c - 

Both theoretical analysis and extensive simulations 
suggest that the signals f(s) which terminate fo~ 
avalanches are uncorrelated and evenly distributed be- 
tween (fo, f c ) provided that f — > f c . The direct conse- 
quence of this observation is that the probability of an 
/o-avalanche merging to fo + ci /o-avalanche is prescribed 
by yz~fo- ^ is important to note that any two subse- 
quent avalanches are mutually independent for the fol- 
lowing arguments to be true. In other words, the proba- 
bility distribution of / -avalanches, initiated immediately 
after the termination of an /o-avalanche of size S is in- 
dependent of S. This is true because in BS model the 
dynamics within an /o-avalanche is completely indepen- 
dent of the particular value of the signals / (s) > f in the 
background that were left by the previous avalanches. 

Here present the master equation. As fo is raised 
by an infinitesimal amount dfo, the probability "flow- 
ing" out of the size distribution of /o-avalanches is given 
by P(S, fo) j^j ) while the probability "flowing" into is 

given by Es^i ^z^PiS - Si Jo). Let fo - f c and 
dfo — > 0, one can write down the master equation as 

(fc-fo)^§^ = -P(S,fo) 
dfo 

S-l 

+ P(Si,fo)P(S-S 1 Jo). (3) 

Si=l 

The first term on the right hand of the equation expresses 
the loss of avalanches of size S due to the merging with 
the subsequent one, while the second one describes the 
gain in P(S, fo) due to merging of avalanches of size Si 
with avalanches of size S — Si. 

In order to investigate the exact master equation it is 
convenient to make some variable changes. Define h = 
— ln(/ c — fo). Therefore, fo = f c corresponds to h = 
+oo. Since in the master equation fo is chosen to be 
close to f c , h varies from a very large number to +oo. 
Due to the variable change the variable h is chosen from 
the distribution P(h) = e~ h , which seems to be more 
" natural" . In the following part we will use the new 
variable h instead of fo- The master equation can be 
rewritten, in terms of h, as 
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dP(S, h) 
dh 



s-i 



= -P(S, h)+J2 p ( s > h ) p ( s - s i> h )- ( 4 ) 



Si=l 



Making Laplace transformation of Eq.(4), after some 
calculation, one obtains 



dHi- P {(i,h)) 

dh 



■ P (f3,h), 



(5) 



where p(f3, h) = J2 P (S, h)e~P s . This exact equation 
s=i 

is the key one of this work. Many interesting physical 
features can be derived from it. As h < +00 avalanches 

size will have a cutoff. The normalization of P(S, h) can 

00 

be expressed as p(0, h) — p (S,h) = 1. Expanding 

s=i 

both sides of Eq. (5) as Tylor series throughout a neigh- 
borhood of the point (5 = 0, one can immediately obtain 

U l - ( S )hP + ^(S 2 ) h P 2 - i(S 3 ) h (3 3 + ...] = 
[(S) h (3 - ±(S 2 ) h f3 2 + ±{S 3 ) h f3 3 + ...] x 
[-l + (S) h p-^(S 2 ) h p 2 + l i (S s ) h ^ + ...] . (6) 

Since the equation (6) holds for arbitrary /?, comparing 
the coefficients of different powers of (3 in the above Tay- 
lor series gives an infinite series of exact equations. Com- 
parison of the coefficients of (3 1 results in 



d\n{S) f 
dh 



= 1. 



(7) 



Eq. (7) is extremely interesting. Changing variable h 
back into /o, one can obtain the "gamma" equation [6,11] 



dfo fc - fa 



(8) 



Inserting the scaling ansatz [6] (S)f ~ (f c — fa) 7 into 
Eq. (8), one immediately obtain an interesting result 



7 = 1 



(9) 



It should be noted that 7 = 1 is universal, that is, 
independent of the dimension. The value of 7 for Jo- 
avalanches is different from those for / -avalanche found 
in Ref. [6], which are 2.70 and 1.70 for ID and 2D 
BS models respectively. Extensive simulations show 
7 = 0.99 ± 0.01 and 7 = 0.98 ± 0.01 for ID and 2D 
BS models respectively. Fig. (1) shows our simulation 
results, which confirms the universal result 7 = 1. 

Higher powers of (i gives new exact equations. Here 
present the first two 



d JS 3 ) h (S 2 ) 



dh { 3(S) h 



2w h ] ^ {s2)h - 



(10) 



(11) 



Next present the solution of the exact gap equation for 
/o-avalanches. Inserting the scaling relation (S)p^ <~ 
(f c — F(s))~ 1 into the equation and integrating, one ob- 
tains 



A/» = f c - F(s) 



where p is the exponent of relaxation to attractor [6]. 
Thus, we obtain p = 1. Interestingly, p is also a universal 
exponent for all dimensional BS models. It shows that 
the critical point (A/ = 0) is approached algebraically 
with an exponent — 1. 

Up to now, we have obtained some exponents of 
corresponding physical properties of /o-avalanches: r, 
avalanche size distribution [10], D, avalanche dimension 
[10], 7, average avalanche size [10], and p, relaxation 
to attractor [6]. Recall another two exponents [6]: v, 
a, which are defined as r co ~ (f c — /o) _ly and 5 co = 
( f c — fo)~v respectively. Here r co and S co are referred to 
as the cut-off of the spatial extent of avalanches (due to 
the limit system size ) and that of the avalanche size (due 
to the fact that / is not chosen as f c ) respectively. It 
is natural to establish some scaling relations of these ex- 
ponents for /o-avalanches similar to those found in Ref. 
[6,12] for /o-avalanches. Nevertheless, these two types 
of avalanches manifest similar fractal properties. Hence 
some common features should be shared by them. In- 
tegrating of the equation (S) = J^~ f ° r ~ SP(S, f )dS 
and the scaling (S) <~ (f c — /o)^ 1 result in 

7= — = 1- (13) 
a 

Due to the compactness [6] of avalanches, we have S co ~ 
rg = (A/)"^, thus 

1 1 



aD (2-t)D' 



(14) 



Eqs. (11)-(12) establish scaling relations among the crit- 
ical exponents, and they imply that the self-organization 
time to reach the critical state is independent of the 
initial configuration of the system. A system of size L 
reaches the stationary state when [A/(s)]" !/ ~ L. It can 
be inferred from Eqs. (11)-(12) that, if one chooses t and 
D as two independent exponents other exponents can be 
expressed in terms of them. Among the six exponents 
mentioned above, r and D can be numerically measured 
[10], 7 and p can be analytically obtained, while v and 
o are difficult to explore despite some methods measur- 
ing the corresponding exponents for / -avalanches were 
introduced in Ref. [13]. Therefore, we can rely on the 
scaling relations and values of the exponents obtained to 
predict the values of v and a. We predict a = 0.2 (ID) 
and 0.275 (2D),_u = 2.04 (ID) and 1.17 (2D). 

Comparing /o-avalanche with /o-avalanche we find 
that the former is more readily to be treated. Two criti- 
cal exponents can be analytically obtained and are found 
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to be universal for all dimensional BS models. Further- 
more, the infinite hierarchy of exact equations and the 
exact gap equations, together with their solutions, pro- 
vide exclusive investigation of the new type of avalanches. 
Another asset of /o-avalanche is that it involves some in- 
formation concerning the whole system. It can be con- 
cluded that /o-avalanche does enable us to comprehend 
the complex system from an effective and different con- 
text. The weak point of this avalanche is that it loses 
some knowledge of individuals. It is still unknown how 
these individual features will matter. It is worthwhile to 
investigate the avalanche dynamics further in the future. 
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Figure Captions 

Fig. 1: The average size of avalanches (S) vs (f c — /o) 
for (a) ID and (b) 2D Bak-Sneppen evolution models. 
The asymptotic slope yields 7 = 0.99 ± 0.01 and 0.98 ± 
0.01 respectively. 
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